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Abstract 

Semi-inclusive deep inelastic scattering of polarized leptons off hadrons enables one to mea- 
sure the antisymmetric part of the hadron tensor. For unpolarized hadrons this piece is odd 
under time reversal. In deep inelastic scattering it shows up as a (sin cj>) asymmetry for the 
produced hadrons. This asymmetry can be expressed as the product of a twist-three "hadron 
— * quark" distribution function and a time reversal odd twist-two "quark — > hadron" frag- 
mentation function. This fragmentation function can only be measured for nonzero transverse 
momenta of the produced hadron. 
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In this paper we study deep inelastic leptoproduction of hadrons. To be precise we investigate 
the process 

£ + H — ► £' + h + X, (1) 

where £ and £' are the incoming and scattered lepton with momenta k and k', H is the hadronic 
target with mass M and momentum P, h is the produced hadron with mass and momentum 
Ph and X is the unobserved rest of the final state. The differential cross section for this process is 
written as a product of a leptonic and hadronic tensor, 

da it a 2 y z 



dx B dzdyd 2 q T 2Q 4 



-L^2MW U , (2) 



where we have used the invariants expressed in the external momenta and the momentum q of 
the exchanged virtual photon (we will limit ourselves to electromagnetic interactions), which is 
spacelike and large (— q 2 = Q 2 3> M 2 ). The invariants used areQ 

Q 2 P-Ph_. 2Ph-q .. Pq (3) 



2Pq P-q Q 2 ' " P ■ k 

The (two-component) vector q T is the transverse momentum of the photon in the frame in which 
the hadrons H and h are collinear. This translates to a perpendicular momentum of the produced 
hadron equal to Ph± = —zQt 1X1 the frame in which the target hadron H and the photon are 
collinear. The leptonic tensor for polarized leptons with helicity A = ±1 is of the form 



L„ u {k, k', A) = [2kX + 2k v k'^ - Q 2 g^ + 2i\e^ pa q p k a ) . (4) 
and the (unpolarized) hadronic tensor is of the form 

2MW,Al,P,Ph) = ^ / (2 ^p0 ( 2 ^(q + P~Px-P h ) 

x(P\J^)\P x ,P h )(P x ,P h \JM\P)i (5) 

where averaging over initial state spins and summation over final state spins is assumed. 

In all of our considerations we will neglect all contributions of order 1/Q 2 , but keep all con- 
tributions of order 1/Q, as it turns out that the time reversal odd fragmentation function, that 
we investigate, contributes at the 0(1/Q) level. We will also not discuss the factorization of the 
process and the inclusion of radiative corrections. We assume that we can use the diagrammatic 
expansion for the hard process as outlined by Ellis, Furmanski and Petronzio (EFP) In this 
approach diagrams are separated in a hard scattering part and soft correlation functions Q. 

The diagrammatic expansion of EFP is useful for deep inelastic scattering, because it in an 
appropriately chosen (color) gauge corresponds to an expansion in powers of Q- 1 . This twist 
expansion as used here is discussed in refs ||, |||, ||, ||. The Born diagrams in this expansion are 
given in Fig. 1 and give || 

2MW flu (q,P,P h ) = ^ J d 4 pd 4 k5\p + q-k)Tr(^(p) lu A(k)) + j Q ^~ <? 1 , (6) 



where the (soft) quark correlation functions 0, ||, || are given by 

1 

(2l0 



*«(P) = ToZaI / d ' xeW ' X (P\^j{0)Mx)\P), (7) 



= £ (2^)4 / d 4 xe^(0\Mx)\Ph,X)(P h ,X\^(0)\0). (8) 



1 With an approximate equal sign we indicate relations valid up to order 1/Q 2 
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In order to analyze the result it is useful to use lightcone coordinates for the external momenta P, 
Ph and q and the integration variables p and k. Using the representation p = \p~ ,p ,Pt\ where 



P 



we write in a frame where H and h are collinear 

' x B M 2 A 



P 



Ay/2 ' x B \f2 
' _ A 

" Q 2 



T 



q T 



A 

AyfT "72 

AV=^ k M 
zQ 2 

Ay/2' zQ 2 V2 



0) 



AMI 



,0 T 



This parametrization is consistent with p + q = k and has been made under the assumption that 
the quark momenta are limited to a hadronic scale of the order of a few hundred MeV. More 
precisely, p 2 , P ■ p and k 2 , Ph - k are of hadronic scale. 

In the deep inelastic limit the process factorizes. For the H — > q part one can integrate over 
p~ . One can then analyze the Dirac content of the projections 0] 



& r \x B ,p T ) = ^ 'dp' Tr($T) 



=x B P+, p T 



(10) 



Using for $ the constraints imposed by hermiticity, parity and time reversal invariance one obtains 
for an unpolarized hadron 



= fi(x B ,p T ), 
M 

— e[x B ,p T ), 



(11) 
(12) 

(13) 



where f\ is the well-known twist-two quark distribution function, in which we have kept the 
p T -dependence and e and f x are the twist-three profile functions discussed in ref. @ and ||, 
respectively. The twist of the profile functions is obtained from the power l/(P + ) l ~ 2 appearing in 
the correlation function. Note that only after integrating over transverse momenta and expansion 
in terms of local operators a comparison with the twist known from the operator product expansion 
can be made. 

Similarly, for the q — > h part one can integrate over k + and analyze the Dirac content of the 
projections 



T) 



1 

Tz 



dk + Tr{AT] 



--P-/z, k 7 



(14) 



As the definition of A involves states \Ph,X), which are out-states, one cannot use constraints 
from time reversal invariance. As a consequence one obtains as the most general Dirac content 



A' 7 1 =D x {z 

^[io- 8 "75] _ 



-zki 



-Zkn 



(15) 
(16) 
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AW = ^E(z,-zk T ), (17) 
A^ = ^DHz,-zk T ), (18) 



A [^ 75 ] = Mh^L H ( Z} - z k T ), (19) 

where D\ is the well-known twist-two quark fragmentation function with the fc^-dependence kept. 
This dependence translates into a dependence on the transverse momenta of produced hadrons. 
Furthermore a time reversal odd twist-two fragmentation function appears, which can only 
be measured when the dependence on the transverse momenta is measured, e.g. in azimuthal 



asymmetries. This function H^- has been discussed in ref. [1C] for polarized leptoproduction as 
a means to probe transversely polarized quarks. At the twist-three level one finds the profile 
functions E (ref. D 1 - (ref. ||) and another time reversal odd function H. 

One can gain some insight in the nature of these time reversal odd functions by considering 
the case that hX would have only one possible channel. Using \Ph,X;out) = e 2lShX \Ph, X;in) 
(see [11]), the constraints from time-reversal then would imply that the functions D\, E and D 1 - 
are proportional to cosJ^x, while the functions and H are proportional to smS^x- 

The twist-two functions are the ones that appear if one evaluates the result for the hadronic 
tensor in leading order. This yields 

2MHV = (Zj?rJ 2^[/lA], (20) 

where P = P — (P • q/q 2 ) q, and P 2 = Q 2 /Ax 2 , while the convolution I[f\ D\] is given by 

I[/lI>l] = J d 2 p ± d 2 k ± 5 2 (p ± + q T -k ± )f 1 (x B ,p ± )D 1 (z,-zk ± ). (21) 
As illustration we will consider the situation in which 

f(x B ,p x ) = f(x B ,0) eM~R 2 HPl), (22) 
D(z, -zk\) = D(z, 0) expi-Rlkl). (23) 

In that case the above integral can be written as 

J2 T}2 t>2 



D i\ = p2 J ff2 ex P 



q T R H R h 
R H + R h 



/i(x B ,0r)£>i(^0 T ). (24) 



Integrating over the perpendicular momenta of the produced hadron, i.e. over q T one obtains the 
familiar factorized form fi(x B ) Di(z), where f\{x B ) = J d 2 px fi(x B ,p T ) is the lightcone momentum 
distribution and D\(z) = / d 2 Ph±Di(z, Phi.) is the fragmentation function. Finally it should be 
noted that everywhere the summation over quark flavors has been suppressed. Together with 
the antiquark part, it can easily be reinstated by adding a summation over flavors of quarks and 
antiquarks weighted with the charge squared (ef) and labeling the functions as and 
respectively. 

The Born contribution to also gives a contribution at subleading order (proportional to 
1/Q). These involve the twist-three functions defined above. At the same order, however, also 
contributions from diagrams with gluons, as shown in Fig. 2 for quarks, must be included. They 
give in addition to the result in Eq. |6| contributions that involve quark - quark - gluon correlation 
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functions. It turns out that they in the gauges A + = (for H — > q part) and A =0 (for q — > h 
part) at the twist-three level contribute through the bilocal matrix elements 

Ftjip) = 7^)4 e^(P\^(0)A T (x)Mx)\P). (25) 

Mtj(k) = -i-^|d 4 x e ik - x {0\AUx)M^)\Ph,X){P h ,X\^(0)\0). (26) 



x 



Using the equations of motion for the quark fields, these can be reexpressed in the quark - quark 
correlation functions $ and A. 

At this point it is appropriate to mention that in order to have a color gauge invariant definition 
one must include a color link in the bilocal quark-quark correlation functions in Eqs |7| and || and 



similarly in the quark-quark-gluon correlation functions in Eq. |25j and 26 . For correlation function 
integrated over k~ and kx one is only sensitive to the nonlocality in the x~ direction. In that case 
a link L(0, x) = V exp[ig Jq ds-A(s)] with a straight path can be inserted in the definition of profile 
functions. Such a link operator becomes unity in the gauge A + = 0. When transverse momentum is 
observed and one only integrates over k~ , one also becomes sensitive to separation in the transverse 
direction, although still x + = 0. In that case one also needs to fix the gauge freedom affecting At- 
This requires a choice of boundary conditions [|l2| for At and correspondingly a particular choice 
of path in the link operator. Although it is not hard to find some way to do this in such a way 
that the link operator reduces to unity after fixing the gauge, it remains to be proven that this 
can be combined with a proof for factorization at the twist three level. 

The full (electromagnetically gauge invariant) result of the diagrammatic expansion up to 
0(1/ Q) reads 



2M W^ u = J d 2 p T d 2 k T 5 2 (p T + q T - k T ) 



Q 



+ y±j^^_tj^El 8xBzf ± Dl 



(P u k±_v + P u k ±u ) /l , 
+l W 8Z M~ h \ XBe ~M h]Hl 



(27) 



where p± is constructed from the vector that in the frame where the hadrons (H and h) are collinear 
is pt = [0, 0,Pt] by subtracting the projection along q, i.e. p± = Pt + (Pt ' Qt / Q 2 ) Q- This vector is 
orthogonal to q. Note, however that neglecting 1/Q 2 contributions, the perpendicular components 
are the same, i.e. p T = p ± . 

In order to investigate the observable consequences of the various terms in the hadronic tensor, 
one needs to consider the contraction with the leptonic tensor. This gives 

v <^IH (29) 

(30) 
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where the azimuthal angles are those between the transverse vectors and the (lepton) scattering 
plane, e.g. p T • q T = \p T \ \q T \ cos<f> p . 

As the most general result one can consider the case where one measures the perpendicular 
direction of the struck quark from the jet direction (p±) and the momentum of one specific hadron 
(h) belonging to this jet, characterized by the lightcone momentum fraction z and the perpendicular 
momentum p^±- One obtains (using k± = — boldmath Ph±/z) 



da (e+H->e'+h+jct+x) 47Ta 2 s 



dx B dydzd 2 Ph± d 2 p± 



Q 4 



y 2 , 

y + 1 



x B fi Di 



+ 2(2 - y) \J\-y cos (f> h 



h±\ 



+ 2(2-y)y/l-y cos, 



zQ 
P±\ 



Q 



y 

x B fi \-D J -- D i 
-^f l D ± + x B [h-x B f L ] Di 



V 

z 



+ 2yy/l - yX smc/) h 



'h± 



2y\J\ -y \ sin^j 



zQ 
P±\ 

Q 



Mx, 



Mx B 



x B e 



m 
M 



x B e 



771 

M 



(31) 



where the H — > q profile functions (fi, f and e) depend on x B and p 2 ± , while the q — > h profile 
functions (D\, D L and Hi~) depend on z and (Ph± — zpj_) 2 . 

The theoretically simplest results are the ones obtained by integrating over the momenta of 
produced hadrons M. Integrating over the perpendicular momenta one has 



da {e+H^e'+h+jct+x) 
dx B dydzd 2 pj_ 



Ana s 



+ 1-2/ 



x B fi(x B ,p ± ) D 1 (z) 



2(2 -y) y/T^ cos <j>j ^ x\f^ (x B , P± )D 1 (z) 



(32) 



This shows that the ratio x B f ± (x B ,p ± )/ fi(x B ,p±) can be obtained most cleanly from the az- 
imuthal (cos (ftj) asymmetry of the produced jet. Note that this ratio is an extension of the naive 
parton result of unity 13]. A calculation in the bag model is shown in ref. [14|. 

The measurement of the time reversal odd fragmentation function is possible from the 
general cross section given above. However, the detection of a jet in the forward direction will be 
hard in view of the fact that one is interested in a twist-three piece (oc 1/Q). This makes it more 
realistic to consider l+H — > £' + h+X. Only convolutions of profile functions from the distribution 
and fragmentation part will appear in this case. We will consider again the approximations for 
the transverse momentum dependence as in Eqs. 22 and and express the results in convolutions 



/ as defined in Eq. 21. The result for the cross section for electroproduction of a hadron h up to 
0(1/Q) becomes 



da ^+H^' + h + X) _ Ana 2 s f 

dx B dydzd 2 p ± Q 4 | 

+ 2(2 - y) \J\-y cos <j> h 
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y 

— + 1 

2 

Ph±\ 



x B I[hDx] 



zQ 



H\ x b 



h(D 1 - -D- 

z 



RhXb ■I[f^D 1 ] 



+ 2y^l - y Asin ( 



h± 



( Mx r 



zQ 



M h R 2 H + R 2 h 



I 



m 
~M 



x B e H 




(33) 
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The time reversal odd fragmentation function is a leading twist fragmentation function which 
appears as a consequence of the fact that the states \ph,X) in the definition of A are out-states. 
This means that it is nonzero because of the interactions of the outgoing hadron h. Furthermore, 
this time reversal odd structure function appears in the cross section as part of a product that 
also contains the profile function e. This chirally odd function has been discussed in ref. § and 
an estimate using the bag model has been given. We note that the combination that appears 
in the cross sections corresponds to a pure quark-quark-gluon correlation function, to be precise 
i(Mx B e-mfi) = (1/2) / dk~ Tr (g F a a a+ ). 

In principle the extraction of the sin<^/j asymmetry for pions does not seem to be extremely 
difficult. One does not need target polarization, while the result is proportional to the lepton 
helicity. The azimuthal structure in the quark transverse momenta plays the key role in this. 
The precise | p± | dependence and dilution of it by QCD corrections (Sudakov effects Jl5|]) are 
less relevant. It is important to emphasize that the pj_-dependence that we have discussed is the 
low-momentum range, say below 1 GeV. Here the p^-dependence is dominated by the intrinsic 
transverse momentum of the profile functions, rather than by perturbative QCD processes [1C] 



that dominate the high transverse momenta 17, 18, 19]. The fact that one deals with a twist three 



structure function furthermore requires not too large values for the momentum transfer Q. It is 
important, however, to have good particle identification and sufficient azimuthal resolution in the 
forward direction. 
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List of Figures 
fig. 1: The Born terms in the expansion of the hadronic tensor. 

fig. 2: The contributions from quark - quark - gluon correlation functions in the expansion of 
the hadronic tensor. Only the quark part is shown. 
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